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ABSTRACT 

Recent calculations by Quataert et al. (2002) found that the growth rates 
of the magnetorotational instability (MRI) in a collisionless plasma can differ 
significantly from those calculated using MHD. This can be important in hot 
accretion flows around compact objects. In this paper we study the transition 
from the collisionless kinetic regime to the collisional MHD regime, mapping out 
the dependence of the MRI growth rate on collisionality. A kinetic closure scheme 
for a magnetized plasma is used that includes the effect of collisions via a BGK 
operator. The transition to MHD occurs as the mean free path becomes short 
compared to the parallel wavelength 2Tr/k\\. In the weak magnetic field regime 
where the Alfven and MRI frequencies u are small compared to the sound wave 
frequency fenco, the dynamics are still effectively collisionless even if uo <C u, 
so long as the collision frequency v <C fciico; for an accretion flow this requires 
v Q\fj3. The low collisionality regime not only modifies the MRI growth rate, 
but also introduces collisionless Landau or Barnes damping of long wavelength 
modes, which may be important for the nonlinear saturation of the MRI. 



1. Introduction 

Balbus & Hawley (1991) showed that the magnetorotational instability (MRI), a local 
instability of differentially rotating magnetized plasmas, is the most efficient source of angular 
momentum transport in many astrophysical accretion flows (see Balbus & Hawley 1998 for 
a review). The MRI may also be important for dynamo generation of galactic and stellar 
magnetic fields. Most studies of the MRI have employed standard MHD equations which 
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are appropriate for collisional, short mean free path plasmas. Recently, however, Quataert, 
Dorland & Hammett (2002; hereafter QDH) studied the MRI in the collisionless regime 
using the kinetic results of Snyder, Hammett & Dorland (1997). They showed that the MRI 
persists as a robust instability in a collisionless plasma, but that at high f3 1 (ratio of 
plasma pressure to magnetic pressure), the physics of the instability is quite different and 
the kinetic growth rates can differ significantly from the MHD growth rates. 

One motivation for studying the MRI in the collisionless regime is to understand radia- 
tively inefficient accretion flows onto compact objects. An example of non-radiative accretion 
is the radio and X-ray source Sagittarius A*, which is thought to be powered by gas accreting 
onto a supermassive black hole at the center of our galaxy (see Quataert 2003 for a review). In 
radiatively inefficient accretion flow models, the accreting gas is a hot, low density, plasma, 
with the proton temperature large compared to the electron temperature (T p ps 10 12 K 
3> T e ps 10 10 — 10 11 K). In order to maintain such a two-temperature configuration, the ac- 
cretion flow must be collisionless in the sense that the timescale for electrons and protons to 
exchange energy by Coulomb collisions is longer than the inflow time of the gas (for models 
of Sagittarius A*, the collision time close to the black hole is ~ 7 orders of magnitude longer 
than the inflow time). 

In this paper we extend the kinetic results of QDH to include collisions; we study 
the behavior of the MRI in the transition from the collisionless regime to the collisional 
MHD regime. Instead of using a more accurate but very complicated Landau or Balescu- 
Lenard collision operator, we use a simpler Bhatnagar-Gross-Krook (BGK) collision operator 
(Bhatnagar, Gross & Krook 1954) that conserves number, momentum and energy. 

There are several reasons for studying the behavior of the MRI with collision frequency: 
(1) one gains additional understanding of the qualitatively different physics in the MHD 
and kinetic regimes, (2) one of the key differences between the MRI in kinetic theory and 
MHD is the anisotropic (with respect to the local magnetic field) pressure response in a col- 
lisionless plasma (QDH). Even if particle collisions are negligible, high frequency waves with 
frequencies ~ the proton cyclotron frequency may tend to isotropize the proton distribution 
function. Our treatment of "collisions" can qualitatively describe this process as well; and 
(3) the transition from the collisional to the kinetic MRI could be dynamically interesting 
if accretion disks undergo transitions from thin disks to hot radiatively inefficient flows (as 
has been proposed to explain, e.g., state changes in X-ray binaries; Esin et al. (1997)). For 
example, there could be associated changes in the rate of angular momentum transport (a). 

The paper is organized as follows. In the next section (§2) we briefly discuss the lin- 
earized kinetic equations in the long- wavelength, low frequency limit (the "MHD" limit); 
this is a review of the formalism used by QDH. In §3 we then derive the kinetic equation 
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for the perturbed pressure including effects of proton-proton collisions via a BGK operator; 
this result is needed to "close" our basic equations and derive the dispersion relation for the 
plasma. In §4 we discuss simpler Landau-fluid (Snyder et al. 1997) closure schemes for de- 
riving the perturbed pressure. The Landau-fluid closure approximations agree well with the 
exact kinetic results from §3 in the low and high collisionality regimes and provide a smooth 
transition for intermediate regimes. In §5 we numerically solve for the growth rate of the 
kinetic MRI and discuss the effects of collisions. Finally in §6 we summarize our results and 
discuss their astrophysical implications. 



2. Linearized Kinetic-MHD Equations 

The analysis in this paper is restricted to fluctuations that have wavelengths much larger 
than proton Larmor radius and frequencies well below the proton cyclotron frequency. In 
this limit, a plasma can be described by the following magneto-fluid equations(Kruskal & 
Oberman 1958; Rosenbluth & Rostoker 1959; Kulsrud 1983): 

^ + V-(pV) = 0, (1) 

p d J + P (V.V)V = ^^-V.P + F g , (2) 

— = Vx(VxB), (3) 
P=P±I+(p\\-p±)bb, (4) 

where p is the mass density, V is the fluid velocity, B is the magnetic field, F g is the 
gravitational force, b = B/|B| is a unit vector in the direction of the magnetic field, and 
I is the unit tensor. In equation (3) an ideal Ohm's law is used, neglecting effects such as 
resistivity. P is the pressure tensor that has different perpendicular (p±) and parallel (py) 
components with respect to the background magnetic field (unlike in MHD, where there is 
only a scalar pressure). The pressures are determined by solving the drift kinetic equation 
given below. P should in general be a sum over all species but in the limit where ion 
dynamics dominate and electrons just provide a neutralizing background, the pressure can 
be interpreted as the ion pressure. This is the case for hot accretion flows in which T p 3> T e . 

We assume that the background (unperturbed) plasma is described by a non-relativistic 
Maxwellian distribution function with equal parallel and perpendicular pressures (tempera- 
tures). Although the equilibrium pressure is assumed to be isotropic, the perturbed pressure 
is not. We take the plasma to be differentially rotating, but otherwise uniform (we neglect 
temperature and density gradients). Equilibrium analysis for equation (2) in presence of a 
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subthermal magnetic field with vertical (B z = B sm9) and azimuthal (B^ = B cos9) com- 
ponents gives a Keplerian rotation (f2 oc R~ 3 ^ 2 ) provided the magnetic field is sufficiently 
weak (B 2 <C GM p /R, where M is the mass of the central object). 

In a differentially rotating plasma, a finite Br leads to a time-dependent B^, which 
greatly complicates the kinetic analysis (unlike in MHD, where a time-dependent B^ can 
be accounted for; Balbus & Hawley 1991); we therefore set Br = 0. For linearization we 
consider fluctuations of the form exp[— iut + ik • x], with k = k R R + k z z, i.e., axisymmetric 
modes; we also restrict our analysis to local perturbations for which |k|i? ^> 1. Writing 
P = Po + 8p, B = B + <5B, p± = po + 5p±, and p\\ = po + Sp\\, V = <fiQR + 5v (with 
Keplerian rotation Q(R)), and working in cylindrical coordinates, the linearized versions of 
equations (l)-(3) become (QDH): 



u5p = p k ■ 5v, (5) 

-iujp 5v R - po^QSv^ = (B Z 5B Z + B^SB^) + %kzBz ^ BR _ ik R 5p±, (6) 

47T 47T 

-iujpoSvj, + p 5v R ^ = l h^_^± _ ifc z s i n g C os 9[5p\\ - Sp±], (7) 

—iujp 5v z = — Lz — - — ik z [sm 2 65p\\ + cos 2 05p±], (8) 

uSB R = -k z B z 5v R , (9) 

uSB^ = -kgBgSvj - lkzBz j^ n Sv R + -B^k • fry, (10) 

ami? 

a;55 2 = k R B z 5v R , (11) 



where k 2 = 4f2 2 + rff2 2 /rflni? is the epicyclic frequency. To complete our system of equa- 
tions and derive the dispersion relation for linear perturbations, we need expressions for 
5p± and 5p». These can be obtained by taking moments of the linearized and Fourier 
transformed drift-kinetic equation that includes a linearized BGK collision operator. The 
drift-kinetic MHD model is described by Kulsrud (1983) based on earlier work by Kruskal 
& Oberman (1958) and Rosenbluth & Rostoker (1959). The drift-kinetic equation for the 
distribution function including the effects of gravity is 

| + (,,b + v.) . V/ + (-6 • - „b • VB + ±(E, + F slM ) I- = C (/) , (12, 

where vg = c(E x B) /B 2 , p = (vj_ — v^) 2 /2i? is the magnetic moment (conserved in our 
approximations in the absence of collisions), F g \\ = GMqwiR ■ h/R 2 , and D/Dt = d/dt + 
^f)|b + v E ^ • V. The fluid velocity V = v E + bu\\, so the Ex B drift v E is determined by the 
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perpendicular component of equation (2). [Only the E x B drift appears directly in equation 
(12). Other drifts such as the grad B, curvature, or gravity xB drifts are higher order in the 
MHD drift kinetic ordering (Kulsrud 1983), which assumes the frequencies are low compared 
to the cyclotron frequency and the gyroradius small compared to gradient scale lengths. On 
the other hand, the parallel component of the gravitational force F g \\ is included as it can be 
the same order as the parallel electric field, which is small compared to the perpendicular 
electric field in ideal MHD.] Note the addition of a collision operator on the right hand side 
of the kinetic equation to allow for generalization to collisional regimes. In the next section 
we derive the linearly-exact kinetic expressions for Sp\\ and 5p±_ using the BGK collision 
operator in equation (12). We then compare these with closure approximations from Snyder 
et al. (1997). 



3. Kinetic Closure Including Collisions 

In this section we use a simple BGK collision operator (Bhatnagar et al. 1954) to cal- 
culate 5p\\ and 5p± from equation (12). Since we consider only ion-ion collisions, the BGK 
operator is Ck (/) = — v (/ — Fm) where v is the ion-ion collision frequency and Fm is a 
shifted Maxwellian with the same density, momentum and energy as / (so that collisions 
conserve number, momentum, and energy). The integro-algebraic BGK operator greatly 
simplifies the calculations while adequately modeling many of the key properties of the full 
integro-differential collision operator. In some situations the effects of weak collisions can 
be enhanced in a more complete collision operator due to sharp velocity gradients in the 
distribution function. We leave investigation of such effects to future work. In this section, 
we calculate the linearization of the drift-kinetic equation around an accretion disk equilib- 
rium including equilibrium flows and gravity. It turns out that a number of complicating 
intermediate terms end up cancelling, and the final forms of the closures used (from equa- 
tions (26-27) onwards) are identical to what one would get from perturbing around a simple 
stationary slab equilibrium. We carried out the more detailed calculation to verify that there 
were no missing terms in the final closures. 

The equilibrium distribution function f is given by 

/o = wikw exp (~^ |v ~ Vo ' 2 ) ' (13) 

where V = v EQ + -Uj| b is equal to the Keplerian rotation velocity in direction. Since 
|v — V | 2 = (v\\ — u\\q) 2 + 2fiB , f can be expressed in terms of as 
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We shall linearize the drift-kinetic equation and the BGK collision operator. The distribution 
function is given as / = /o + Sf where Sf is the perturbation in the distribution function. 
The shifted Maxwellian that appears in the BGK collision operator is given by 



M 



F M has three free parameters (N M , u\\m, T m ) which are to be chosen so as to conserve 
number, parallel momentum and energy. When taking moments of the BGK operator, it is 
important to note that J d 3 v = f 2ir (B + SB) d/idv\\. From equation (15) and conservation 
of number, momentum and energy it follows that 



N, 



= n + Sn ps n ^1 + + 2^.00 J d/idv\\5f, (16) 

NmU\\m = N M (u\\ + Su) ps n u\\ ^1 + + 2ttB J d/idv^Sfv^ (17) 
N M T M = Po + 6p = p + (S Pll +2Sp ± )/3, (18) 
5p\\ ~ Po5B/B + 2nB J dfj,dv\\5fm(v\\ — u\\ ) 2 , (19) 

Sp± ~ 2p 5B/B + 2nB J dfxdv l{ 5ffimB , (20) 

where the approximate expressions retain only linear terms in perturbed quantities. Lin- 
earizing the expression for the relaxed Maxwellian in equation (15) about f , the drift-kinetic 
BGK collision operator is given by 

C K (Sf) = -uSf + uf x 

({Sn 3ST\ m(, w , \i ST\ mfiB fSB ST\).^. 

The drift-kinetic equation including the BGK operator can be linearized to obtain the fol- 
lowing equation for Sf 

., . w . .(SB^smO - 5B Z cos 9) m m (y\\ - u\\ ) f 

Sf = U4o(y\\ -u ll0 )f sm9 ? — + — , \" X 

io-oo io + ik\\ [v\\ - u\\ ) + v) 

+ M±M) + 7 . ^ . . x 

m I [—luj + ik\\ [v\\ — u\\o) + v) 

Sn 3 5T ^ m(v\\ - u\\ )5u ^ m{v\\ - up) 2 ST ^ m(j,B ST m(j,6B\ 
n 2 T T 2T T T T T / ' 

where F s y = GM mSB R /B R 2 is the component of gravitational force in the direction of 
magnetic field. Choosing a compact notation where — icu sin (&B<£ sin 9 — SB Z cos #) mu^ / eB 
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+F g \\/e + E\\ — > E\\, the moments of the perturbed distribution function 5f in drift coordi- 
nates (v\\,fj), f (l,2p,B , (v\\ - up) 2 ) Sf2nB d[idv\\ give 



n B y ' ik\\T \\n Q 2 T J \T B J 



(23) 



Po B y ' ik\\T {\n 2T J \T B J 

+V2^R + + 2jsin6 ,V^(^sing-^cosg)\ j 
Co V^o ^ ^0 /J 

+ V2- (1 + 2 C 2 i?) + f ^ + 2, sin /ll^o(^^-^cos0) \ ( + 1 

Co V-M) f B J J 

£"11 can be eliminated by taking appropriate combinations of these three equations: 

*> ^ = J4 {1 - R)+C2Z (^- S -i-), ,26) 



Po Po B \T B J 
and 

(1 + 2C*) £ - «^ = f (1 + 2C 2 fl - fl) - 6 - - f - f ) , (27) 

Po Po -t>o \Po ^0 -DO/ 

where <5T = (25T X + 5T N ) /3, 55 = b • 5B, C = + w) /y/2\k\\\co, (2 = iv/y/2\k\\\co, 
fey = b • k, T|| X = mp\\ t ±/ p, and c = ^/T /m is the isothermal sound speed of the ions. In 
equations (26) and (27), R = 1 + is the plasma response function, where 

^=T«h^ (28) 

is the plasma dispersion function (NRL plasma formulary 2000). Equations (26) and (27) 
can be substituted into the linearized fluid equations in §2 to derive the dispersion relation 
for the plasma. The full closures are, however, very complicated, so it is useful to consider 
several simplifying limits that isolate much of the relevant physics. In addition, the solution 
of the MHD equations from §2 with fully kinetic closures will give an implicit equation for 
the growth rate (involving the Z function) that would have to be solved numerically. 

The closure equations can be simplified in two limits, |£| <C 1, the collisionless limit, 
and \(\ 3> 1, the high collisionality limit. The derivation of the asymptotic solution for the 
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closure equations in these two limits is given in the Appendix. For high collisionality 

sp, sip + <i(i + syj>_ 2 w (29) 



Po 3p C2 V 3 9 CfV Po C2 B ' 
and 

*!L = 5 S _R , Ci f_2 5 \ 5_p Ci 55 

Po 3 po C2 V 3 + 9dV Po C2 A) ' 1 J 

where Ci = ^/V^\k\\\co- Notice that in the limit that the collision frequency is very high, 
(2 — > 00, and one recovers the MHD result that the perturbations are adiabatic and isotropic: 
fip\\/Po = $P±/Po — 55p/3p . Inspection of equations (29) and (30) suggests that the MHD 
limit will be reached whenever IC1/C2I <C 1, i- e -j v ^ Later we shall show that in fact 
v ^> v^2 1 A: 1 1 1 Co is required, i.e., the collision time must be much less than the sound crossing 
time of the wavelength of the mode. This is important because the MRI has u <C |&|||co in 
a high (5 plasma so the regime uj <C v <C \k\\\co is an interesting one. 

For low collisionality, \(\ <C 1, to second order in (, 

S P± S P ■ 1 — /■ SB 7r CiC2 Sp t t ( n\SB 

+C1C2 2-- — , 31) 

Po 

and 



Po Po B 3 po V 3/ B t 



~ WTrCl -5-)+—\ 4 ClC2 - TtCi 7, + 



Po Po VPo Bo J Po 



Wq [yfiUM* - - 2C 2 - 4C2CJ • (32) 

To first order, there is no effect of collisions on the growth rate of the MRI; the results above 
are then exactly same as equations (20) and (21) in QDH (who neglected collisions entirely). 
Collisional effects modify the closure only at order ( 2 , though one has to go to this order to 
find the first order dependence of uj on v in the dispersion relation. 



4. Comparison with Landau-Fluid Closures 

The results from the last section provide useful expressions for 5p± and 5p\\ in the low 
and high collisionality regimes, \(\ <C 1 and \(\ 3> 1, but it would be convenient to have a 
single set of equations that can provide a robust transition between these two regimes. The 
Snyder et al. (1997) closure approximations, which we discuss in this section, can do this. 

The second moments of the drift kinetic equation (eq. [12]) yield evolution equations 
for 5p± and 8p\\ (see, e.g., eqs. [16]- [17] of Snyder et al. 1997). The linearized versions of 
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these equations, including a BGK collision operator, are given by 1 

— iu8p\\ + poik. ■ 5\ + ik\\q\\ + 2pQik\\8v\\ — 3po^ cos — = —-v {bp\\ — 5p_i) , (33) 

i?o 3 

and 

— iuj5p± + 2p i\i ■ 5v + ih\\q± — p ik\\8v\\ + -Po& cos ^ R = — v [8p± — 8p\\) . (34) 

2 Bq 3 

As is usual with moment hierarchies, the above equations for 8p\\,± depend on third moments 
of the distribution function, q\\ and q±, the parallel and perpendicular heat fluxes. Snyder 
et al. introduced closure approximations for q» and q± that determine 5p± and 8p» without 
solving the full kinetic equations of the previous section. These Landau-fluid approximations 
"close" equations (l)-(4) and allow one to solve for the linear response of the plasma. 

The linearized heat fluxes in the perpendicular and parallel directions are given by 

2 ik\\(Sp ± /po-Sp/p ) 

q± = ~PoC , — r- (35) 

(^V7r/2|A;|||co + v] 

and 

n = ^ A^^ - . (36) 

y " y (V8^\k l{ \c + (3tt- 8) v) 1 ; 

As discussed in earlier work (Snyder et al. 1997; Hammett et al. 1992, 1993; Smith 
1997), Landau-fluid closure approximations provide n-pole Pade approximations to the exact 
plasma dispersion function Z(() that appears in the kinetic plasma response of §3. These 
Pade approximations are thus able to provide robust results that capture kinetic effects 
such as Landau damping, and that can also smoothly transition between the high and low 
( regimes. 2 We have found that, not surprisingly, the fluid approximations remain robust 
when collisions are included. That is, in all of the numerical tests we have carried out, we 
have found good agreement between the results from equations (33)-(36) and the asymptotic 
kinetic results from the previous section for the low and high collisionality regimes. Thus all 



1 A comparison of our equations (33) and (34) with equations (30) and (31) in Snyder et al. shows that 
our equations have an extra term proportional to the Kcplcrian rotation frequency; this is because Snyder 
et al. (1997) did not include gravitational effects and Keplerian rotation in their linearized equations. 

2 Thc approximations are fairly good near or above the real ( axis, though they will have only a finite 
number of damped roots, corresponding to the finite number of poles in the lower half of the complex plane, 
while the full transcendental Z(() function has an infinite number of damped roots. 
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of the plots in this paper are calculated with the Snyder et al. (1997) Landau-fluid closure 
approximations of equations (33)- (36). 

The Snyder et al. Landau-Fluid closure approximations provide a useful way to extend 
existing non-linear MHD codes to study key kinetic effects. The closure approximations are 
independent of frequency (or the Z function), and so are straightforward to implement in an 
initial value nonlinear code (though they do require FFT's or non-local heat flux integrals to 
evaluate some terms(Snyder et al. 1997; Hammett et al. 1992)). But one should remember 
that they are approximations and so do not accurately model all kinetic effects in all regimes, 
particularly near marginal stability (Mattor (1992); Smith (1997); Dimits et al. (2000)), 
though we have generally found in other applications that they work fairly well in strong 
turbulence regimes (Hammett et al. (1993); Parker et al. (1994); Smith (1997); Dimits et 
al. (2000)). 

As an aside, we note that the double adiabatic theory of Chew, Goldberger, & Low 
(1956), which is a simpler closure approximation that sets q\\ = q± = in equations (33) and 
(34), generally does a poor job of reproducing the full kinetic calculations. This is because 
the perturbations of interest have uj <C | fcy |c and are thus far from adiabatic (see also QDH). 

5. Collisionality dependence of the MRI growth rate 

Figures 1 and 2 show the growth rate of the MRI for intermediate values of collisionality 
in addition to the limits of zero and infinite (MHD) collision frequency (the latter two cases 
were shown in QDH). To produce these plots, we have used equations (5)-(ll) and (33)-(36). 
These equations were solved both with a linear initial value code to find the fastest growing 
eigenmode, and with Mathematica to find the complete set of eigenvalues uj. 

Figures 1 and 2 show that the transition from the MHD to the collisionless regime is 
fairly smooth and occurs, for these particular parameters, in the vicinity of v/VL ~ 10 3 , which 
corresponds to v ~ 10kc , or k\ m f p ~ 0.1, where \ m f p = c /u is the mean free path. Figure 
3 shows the growth rate versus collisionality for f3 z = 100 and f3 z = 10 4 , and for = B z , 
k R = and = 0, k R /k z = 0.5. 

It is clear from these figures that the transition from the collisionless to the collisional 
MRI takes place at far higher collision rates than v ~ f2 ~ uj. That is, v > uj is not a 
sufficient criterion to be in the collisional regime. Instead, the collisional regime requires 
v ^> |A;|||co, which can be written as u/Q ^> \^j3\k\\ \va z /Q ~ y^j3. Figure 3 shows that the 
much of the dependence on collisionality for different values of f3 z can be captured by plotting 
the growth rates vs. v j ' y/~[J~ z VL, though there is some residual variation. 
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At high > 1, the Alfven and MRI frequencies are small compared to the sound 
wave frequency, and there exists a regime uj <C v ^ k\\c where the collisionless results still 
hold despite the fact that the collision time is shorter than the growth rate of the mode. 
Physically, this is because in order to wipe out the pressure anisotropy that is crucial to the 
MRI in a collisionless plasma (see QDH) the collision frequency must be greater than the 
sound wave frequency, rather than the (much slower) growth rate of the mode. This can 
also be seen by comparing Figures 1 and 2 with the corresponding figures in QDH: the effect 
of increasing collisions (decreasing pressure anisotropy) is similar to that of decreasing (5 Z 
(decreasing pressure force relative to magnetic forces). From the point of view of Snyder 
et al.'s fluid approach, the weak dependence of growth rate on collisionality even if v is 
as large as uj can be understood from the fact that the terms proportional to uo and v in 
equations (33) and (34) are both much smaller than the dominant terms involving convection, 
heat conduction, and magnetic forces. So the relative magnitudes of uj and v are not that 
important, and it is not until v is large enough to be relevant in equations (35)-(36) that 
collisional effects become noticeable. 

Figure 4 shows the complete spectrum of eigenmode frequencies as k z is varied, including 
the propagating and damped modes in addition to the unstable MRI branch. We show all 
of the waves present in collisionless Landau fluid and MHD calculations for a fairly general 
choice of wavenumbers and a moderate (3 Z = 10. The MRI is operational at lower k Z) while 
at high k z the eigenfrequencies eventually approach the uniform plasma limit. 

Focusing first on the MHD solutions at high k z , we see the standard set of 3 MHD waves: 
in order of descending frequency these are the fast magnetosonic wave, the shear Alfven 
wave, and the slow wave. Equations (5)-(ll) with an MHD adiabatic pressure equation 
uj5p = p k ■ 5v is a set of 8 equations with 8 eigenvalues for uj. The standard 3 MHD waves 
provide 6 of the eigenvalues (±uu for oppositely propagating waves). The remaining roots are 
zero frequency modes (not shown in the plot). One is an entropy mode, corresponding to 
Sp/po = —5T/T so that Sp = 0. The other solution corresponds to an unphysical fluctuation 
that violates V • B = 0, which is eliminated by imposing the proper initial condition that 
V • B = 0. At lower k z in the MHD plots in Figure 4, it is the slow mode that is destabilized 
to become the MRI, as discussed in Balbus & Hawley (1998). 

Turning next to the collisionless limit in Figure 4, there are two roots plotted in addition 
to the three "MHD-like" modes; this is because the single pressure equation of MHD is 
replaced by separate equations for the parallel and perpendicular pressure, so that there are 
now two entropy-like modes, both of which have non-zero frequencies but which are also 
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strongly damped by collisionless heat conduction (which is neglected in MHD). 3 

The fast, Alfven, and slow waves in the collisionless calculation can again be identified in 
order of descending (real) frequency at high k z . At lower k z , one of the slow modes becomes 
destabilized to become the MRI, as in MHD. Unlike in MHD, however, the fast magnetosonic 
waves are strongly Landau damped since the resonance condition uj ~ k\\Co is easily satisfied. 
In addition, it is interesting to note that both the shear Alfven and slow waves have some 
collisionless damping at the highest k z used in this plot, though the damping will approach 
zero for very high k z . In a uniform plasma the shear Alfven wave is undamped unless its 
wavelength is comparable to the proton Larmor radius or its frequency is comparable to the 
proton cyclotron frequency (neither of which is true for the modes considered here). By 
contrast, the slow mode is strongly damped unless k± <C k\\ (e.g., Barnes 1966; Foote & 
Kulsrud 1979). The damping of small k z shear Alfven waves in Figure 4 is due to the fact 
that our background plasma is rotating so the uniform-plasma modes are mixed together. 
Thus the well-known dissipation of the slow mode by transit-time damping also leads to 
damping of what we identify as the shear Alfven wave (based on its high k z properties). 



6. Summary and Discussion 

In this paper we have extended the linear axisymmetric kinetic magnetorotational in- 
stability (MRI) calculation of QDH to include the effect of collisions. The MHD limit is 
recovered when the mean free path is short compared to the MRI wavelength, i.e., v ^> k\\CQ, 
with a fairly smooth transition between the collisionless and collisional regimes. Interest- 
ingly the collisionless MRI results hold not only if v <C uj, but even when uj <C v <C k\\CQ. 
This intermediate regime can exist in /3 ^ 1 plasmas because the MRI growth rate is slow 
compared to the sound wave frequency, uj ~ k\\VA = k\\c$yj2/ j3 <C k\\c . 

If we consider the application of our results to accretion flows, the collisionless limit will 
be applicable so long as is/Q & y/]3. This condition is amply satisfied for proton-proton and 
proton-electron collisions in all hot radiatively inefficient accretion flow models, suggesting 
that the collisionless limit is always appropriate. However, high frequency waves such as 



3 We should point out that while our equations using the Snyder et al. (1997) 3+1 Landau-fluid closure 
approximations have 8 eigenfrequencies, the equations using the more accurate 4+2 Landau-fluid closure 
approximations have 10 eigenfrequencies, with 2 additional strongly damped roots. If the exact kinetic 
response were used, one would find an infinite number of strongly damped eigenmodes because the Z(() 
function is transcendental. These strongly damped modes are related to "ballistic modes" and transients in 
the standard analysis of Landau damping. 
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ion-cyclotron waves can isotropize the proton distribution function and thus provide an 
effective "collision" term crudely analogous to that considered here. It is difficult to estimate 
the importance of this process (e.g., whether its effective collision frequency is <; Q^/~/3) 
because we don't know to what extent high frequency waves will be excited in the accretion 
flow. They are probably not significantly excited by the underlying MHD turbulence that 
drives accretion since this maintains low frequencies throughout the turbulent cascade (see 
Quataert's 1998 discussion based on Goldreich & Sridhar 1995). High frequency waves may, 
however, be excited by shocks, reconnection events, or velocity space anisotropies. 

One might anticipate that the linear differences between the collisionless and collisional 
MRI highlighted here and in QDH will imply differences in the nonlinear turbulent state in 
hot accretion flows (see, e.g., Hawley & Balbus 2002; Igumenshchev et al. 2003 for MHD 
simulations of such flows). Not only are there differences in the linear growth rates of the 
instability that drives turbulence, but the spectrum of damped modes is also very different. 
In particular, in the kinetic regime there exist modes at all scales in |k| that are subject to 
Landau/Barnes collisionless damping, while in the MHD regime the only sink for turbulent 
energy is due to viscosity/resistivity at very small scales (very high |k|). Indeed, as we have 
shown, even long wavelength Alfven waves can be damped by collisionless effects because 
of the mixture of uniform-plasma modes in the differentially rotating accretion flow (§5 and 
Fig. 4). Whether these differences are important or not may depend on how efficiently 
nonlinearities couple energy into the damped modes. These could modify the nonlinear 
saturated turbulent spectrum (e.g., the efficiency of angular momentum transport) or the 
fraction of electron vs. ion heating (the heating may also be anisotropic), which in turn 
determine the basic observational signatures of hot accretion flows (the accretion rate and 
the radiative efficiency). One approach for investigating nonlinear collisionless effects would 
be to extend existing MHD codes to include anisotropic pressure, the Snyder et al. (1997) 
fluid closure approximations for kinetic effects, and the BGK collision operator considered 
here. By varying the collision frequency, one can then scan from the collisionless kinetic to 
the collisional MHD regime, and assess any differences in the nonlinear turbulent state. 

We thank Drs. W. Dorland, P. B. Snyder and W. Tang for useful discussions. This 
research was supported by NASA grant NAG5-12043, D.O.E. Contract No. DE-AC02- 
76CH03073, NSF grant AST-0206006, and an Alfred P. Sloan Foundation Fellowship (to 
EQ). 
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A. Closure for high collisionality: \(\ >> 1 

For |C| » 1, Z (C) « -1/C - 1/2C 3 - 3/4C 5 , R » -l/2( 2 - 3/4C 4 , 1 + 2( 2 R » -3/2C 2 - 
15/4C 4 , Z - 2(R ss 1/C 3 + 3/C 5 . Equation (26) then becomes 

Sn S P± SB (x 1 x_ ( / 1 1 )( st «y (M) 



no Po B V 2CV C V 2CV V^o B J 

Assuming IC1/C2I <S 1 (a high collisionality limit 00 v) and using the binomial expansion 
we get 



"0 Po I C2 CI \2 "V C 2 V 2 7 J VC2A) T 

To the lowest nonvanishing order one gets 

Sn(i 6p ± (I + 2 Ci\ + 6 _n(\_Ci_\ = _Ci -jB (A3) 



^0(2 Po \3 3(2/ Po \3 3( 2 J (2 B 
Expanding equation (27) gives 

n V 2C 2 + 4CV + po V 2 C 2 4CV ~ fl U 2 C 4 
/5£ 5n 5T V / 1 3 \ 

H*-* + ^M? + cO- (A4) 



This preprint was prepared with the AAS IATgX macros v5.0. 
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Again using the binomial expansion for IC1/C2I *C 1 we get 

5n ( 3Ci , 9 /CiV , 3 \ 5pn A 1 Ci 1 /Ci X ' 



n l 2C2 + 2 VC 2 y + 4C|J + Po I 3 2C 2 + 2VC 2 ; + 4C| 



The lowest order solution is 

2( 2 n \3 2( 2 ) Po V 3 ( 2 ) Po C2 A)' { ' 

We shall expand the parallel and perpendicular pressure perturbations as 5p± = S°p± + 
Ci/C25V+(Ci/C2) 2 5V+- and<5p|| = 5 p\\+(i/(2S 1 P\\ + ((i/(2) 2 S 2 P\\+- From equations (A2) 
and (A5) one gets S°p\\/p = S°p±/p = 5Sn/3n for the lowest order, and (Sp± — S 1 p\\)/p = 
35B/B — 25n/n . To the next order we can expand the solution as 

Sp\\ _ 55n | Ci^P|| | / CiV Pp\\ (A7) 



Po 3n C2 Po VC2/ Po 

5P± Mn + ^ f*lE!L + 3^-2^Vf-V— ■ (A8) 



Po 3n C2 V Po 5 7io/ VC2/ Po 
To the next order in C1/C2 m equation (A2) one gets 

2^ n 2 po 3 V Po Po J 1 ' 

To the next order in equation (A5) we get 

2 >ft.^) ». (A10) 

3Ci / ™o Po 3 V Po Po / B 
Equations (29) and (30) follow from equations (A9) and (A10). 



B. Closure for low collisionality: |£| •< 1 

This regime is useful for low collisionality v <C k\\Co and high (3, where the MRI is low 
frequency as compared to the sound wave frequency Using the asymptotic expansion for 
ICI < 1, Z (C) « - C 2 ) - 2C and R(() ta 1 + iy/n - 2( 2 , we simplify equation (26) to 

get 

5n 5p ± SB , r- . (ST SB\ > , ^ , x 
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The lowest order term in ( gives Spj_/po = 5n/n . Let Spj_/po ~ 5n/n + (5 1 p±/po. To the 
next order one gets 

C— = -»v^C-5- + ^2-5- = -»v^Ci-5-- (B2) 

£>o -do -do -do 

Therefore to second order in (, Sp±/p ~ 5n/n — iy/ir^SB / B + ( 2 5 2 p±/p . On using the 
asymptotic formula for Z and R in equation (27), one gets 

Sn , Sp\\ /—^SB _ , r- k , f 8n 5T» SB\ , . 

— - 1 + iV^C) = -isftC W ~ C2 (iyft ~ 4C — - - -jr ■ B3 
n Po -Bo \ n o 2T j3 / 

To the lowest order one gets 8p\\/po = Sn/no, so let Sp\\/po ~ Sn/rio + C^Pw/Po- To the next 
order, 

C— ^ = -iy/^Ci— + B4 

Po ™0 -»0 

Therefore through second order 5p\\/po ~ Su/uq + iy^Ci (8B/B — 5n/n ) + C 2 ^ 2 P||/Po- The 
comparison of the terms of the order £ 2 in equation (26) give 



= 2C1C v - ttCiC 2 v + - ' (B5) 

Po B 3 V^o no J 

and the terms of the order ( 2 in equation (27) give 

C 2 — = ( 4 Cl ( 2 - <rf - £ d&) - + UdC - - 2C 2 - 4«) f ■ (B6) 

Po V 6 / n \ 6 / j3 

From equations (B5) and (B6) the asymptotic expansion in equations (31) and (32) follow. 
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Fig. 1. — Growth rates of the MRI as a function of kn/k z for different collision frequencies; 
(3 Z = 10 4 and (a) = 0, (b) = B z . For u/Vl > 10 4 the growth rates are very close 
to the MHD values, while for v/Q < 10 2 they are quite similar to the collisionless limit. 
The enhancement of the growth rate in the collisionless regime for small k R is the result of 
pressure anisotropy. 
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Fig. 2. — Growth rates of the MRI as a function of k z for different collisionalities. 
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Fig. 3. — Variation of the MRI growth rate with collisionality for k R = 0, = B z (top 
curves) and k R /k z = 0.5, = (bottom curves). Collisions isotropize the distribution 
function and can increase the growth rate in some regimes and decrease it in others. Con- 
tinuous lines correspond to (3 Z = 10 4 and dotted lines to [3 Z = 10 2 . The transition from 
collisionless to MHD regime takes place when i//Q ks \/]3. 
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damped slow: d, growing slow: g, entropy modes: e1 ,e2 



Alfven modes: A1 ,A2, fast modes: f1 ,f2 




Fig. 4. — The real and imaginary parts of the mode frequency as a function of k z using 
collisionless Landau fluid closures (a,b) and MHD (c,d) are shown [y = 0, k R VA z /Q = 0.5, 
(3 Z = 10, = 0). 



